1 



TIMING WITH THE EPIC PN CAMERA OF XMM-NEWTON 
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Abstract 

The EPIC pn CCD camera on board of XMM-Newton 
is designed to perform high throughput imaging and spec- 
troscopy as well as high resolution timing observations 
in the energy range of 0.1-15 keV. A temporal resolution 
of milliseconds or /i-seconds, depending on the instru- 
ment mode and detector, is outstanding for CCD based 
X-ray cameras. In order to calibrate the different observ- 
ing modes of the EPIC pn CCD, XMM-Newton observa- 
tions of the pulsars PSR B1509-58, PSR B0540-69 and 
the Crab were performed during the calibration and per- 
formance verification phase. To determine the accuracy of 
the on board clock against Coordinated Universal Time 
(UTC), PSR B1509-58 was observed simultaneously with 
XMM-Newton and RXTE in addition. The paper summa- 
rizes the current status of the clock calibration. 
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1. Introduction 

The EPIC instruments on board of XMM-Newton^ suc- 
cessfully launched on 1999 December 10, are designed for 
imaging, high throughput spectroscopy, and timing anal- 
ysis. The EPIC pn camera offers six different observa- 
tion (readout) modes, four imaging modes, and two fast 
readout modes. The imaging modes are dedicated to high 
throughput imaging and spectroscopy, while the fast read- 
out modes are designed for high temporal resolution down 
to 7.2 /is in combination with high spectral resolution. 
For a detailed description of the EPIC pn readout modes, 



their implementation, and characteristics see e.g. Kend 



ziorra et al. (1997), Kendziorra et al. (1999), Kuster et al. 



(1999), or Ehle et al. (2001) 



This paper is structured as follows: In section^ we give 
a brief description of the fast readout modes available for 
the EPIC pn camera and their technical implementation 
as far as these are of importance for the observer. In sec- 
tion 1^ and ^ we present preliminary results of the relative 
time calibration based on observations of young millisec- 
ond pulsars. In section || we describe limiting constraints 



an observer has to take into account when using the EPIC 
pn camera for observations with high temporal resolution. 

2. Fast readout modes 

In addition to the imaging modes the design of the EPIC 
pn camera offers two observational modes designated to 
temporal analysis of X-ray sources with high time resolu- 
tion, called Timing and Burst Mode. While Burst Mode is 
designed for very bright sources up to 6.3 Crab only, Tim- 
ing Mode is usable for faint sources as well. In both modes 
spacial information in y-direction is lost due to the con- 
tinuous readout of the CCD. In contrast to Timing Mode 
which uses the full size of the CCD, the last 20 lines of 
the CCD are discarded in Burst Mode. To improve time 
resolution to the maximum possible, the position of the 
source on the CCD has to be taken into account. This 
is done within the XMM Science Analysis System (SAS) 
for the nominal source position by default. For pointings 
deviating from the nominal position, the observer has to 
specify the source position on the CCD during the extrac- 
tion process. 

Beside high temporal resolution, both modes offer the 
full spectral resolution of the EPIC pn camera. For a re- 



port on t he status of energy calibration see [Briel et al 
1(2002)1 or pKirsch et al. (2002)| in this volume. 



3. Timing accuracy 

All XMM-Newton data is tagged with a time stamp (on 
board time OBT) from a 1 Hz on board clock provided 
by a temperature compensated oscillator of the XMM- 
Newton Command and Data Management Unit (CDMU) 
(see Aranci et al. 1998| ). The EPIC pn event analyzer 
EPEA further uses this clock pulse as an input for its 



Table 1. Mode specific parameters for the EPIC pn camera. 
For the maximum possible time resolution the position of 
the source and the PSF has to be taken into account. 



Obs. Mode 


Frame 


Life 


Time 




Time [ms] 


Time [ms] 


resolution 


Small Win. 


5.672 


4.028 


5.67 ms 


Timing 


5.965 


5.912 


29.56 /is 


Burst 


4.345 


0.126 


7.2 /is 
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Table 2. Observations used for time calibration. 



Object 


Rev. 


Mode 


Chopper 


Obs. Time 

[ksec] 


Crab 


0056 


Timing 


25 


7.0 




0234 


Burst 





10.0 


B1509-58 


0137 


Timing 





9.6 




0137 


Small Win. 





11.4 


B0540-69 


0085 


Timing 





7.0 




0085 


Small Win. 





7.0 



1.0 r 



internal clock which has an accuracy of 16 jis. All events 
handled by the EPEA are tagged with this accuracy. A 
verification of the internal timing of the EPEA was done 
during ground calibration at th e PANTER facility before 
launch ( Kendziorra et al. 1997 ). 

On board time (OBT) is further converted to Coordi- 
nated Universal Time (UTC) on ground while extracting 
data using information stored in a time correlation file, 
delivered with each Observation Data File (ODF). This 
file contains OBT versus UTC information. To transform 
event times, OBT is correlated versus UTC using a poly- 
nomial fit. The user can influence the accuracy of the fit- 



ting routine; for details see Lammers (2001) 



During the extraction process with the SAS, the event 
times are corrected according to the readout mode used 
during observation. In imaging modes the photon arrival 
times are transformed to the center of integration interval, 
for the best time resolution possible. In the fast readout 
modes the time resolution can be improved by taking the 
position of the source and the line in which the photon 
was detected into account. By default line 190 is used as 
source position, which is equal to the nominal pointing 
position. 

The last correction that has to be applied is a trans- 
formation to the solar system bary centre. This is done via 
the SAS tool "barycen" which transforms the given event 
times stored in the event file to Barycentric Dynamical 
Time (TDB). 

4. Observations and data analysis 

During commissioning and performance verification phase 
of X MM-Newton observations of voung millisecond pul- 
sars wi th the purpose of relative and absolute clock cal- 
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Figure 1. Pulse profiles of the Crab pulsar using RXTE 
PCA and XMM-Newtons^ EPIC pn camera in Timing and 
Burst Mode. For RXTE PCA the pulse profile in the en- 
ergy range 2.0-40 keV is shown. For each pulse profile the 
maximum count rate is normalized to unity. Pulse phase 
zero is defined by the maximum count rate in each pulse, 
therefore absolute phase information is lost. 



background to a minimum we ignored flux below 3keV. 
This is of importance especially for observations during 
early commissioning phase, e.g. the Crab observation in 
Rev. 56, when the detector setup was not yet optimized. 
For observations in fast modes we used line 190 as posi- 
tion of the source in RAWY direction for time correction. 
In addition we selected 4-5 columns around maximum in- 
tensity in RAWX direction as extraction region, to reduce 
non pulsed flux from nebula emission or background. Fur- 
ther we transformed photon arrival times to TDB, using 
the SAS tool "barycen" and object coordinates given in 
Tab. |. 

From these light curves we derived pulse periods for 
each object and observation mode using a -maximization 
test. The resulting pulse periods Pfound are given in Tabjl 



Using Jodrell Banks monthly data base for the Crab (Lyne 
et al. 2001) and Princeton Pulsar Data base ( p?aylor et al 



ibratio h were made. In addition we scheduled a simul- 
taneous observation of PSR B1509-58 with RXTE with 
purpose of absolute time calibration. A summary of all 
observational data used for the analysis presented in this 
paper is given in Tab. ^ For all Crab observations in Tim- 
ing Mode the electronic chopper was set to 25 to reduce 
telemetry rate. All other observations were performed with 
electronic chopper set to 0. Using the standard SAS proce- 
dure (SAS version 5.1.0) we extracted event files and light 
curves for all observations, objects and detector modes 
in the energy band 3.0-12.0 keV. To reduce low energy 



1993), we extrapolated radio ephemeris Pradio for all sources 
using the center of observation time of the X-ray observa- 
tions as reference point. Subsequently we folded the X-ray 
light curves with the respective pulse period to get pulse 
profiles shown in Fig. |] and Fig. H. For all pulse profiles 
the maximum flux is set to unity and pulse phase zero is 
defined by the maximum flux. Therefore Fig. |^ and Fig. || 
do not contain absolute phase information. 
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Figure 2. Pulse profiles of PSR B 1509-58. The profiles are 
treated in the same way as the pulse profiles shown in 
Fig. 0. Note that the small shift; in phase between Tim- 
ing and Burst Mode data is not a difference in absolute 
phase. 



5. Observational results 

In Fig. |l| the resulting pulse profile for the Crab is shown 
for Timing Mode, Burst Mode and RXTE PCA obser- 
vations. The results for PSR B1509-58 in Timing Mode, 
Small Window Mode and RXTE PCA are shown in Fig. ||. 

As Fig. and Fig. || clearly demonstrate, we get con- 
sistent pulse shape for all observational modes of XMM- 
Newton. In addition the pulse profiles are well resolved 
and are in agreement with those we get from RXTE PCA 
data. As described earlier we took advantage of the imag- 
ing capability of XMM-Newton and selected only the cen- 
tral emission region in RAWX direction during data ex- 
traction in the fast readout modes. This reduces non pulsed 
flux from nebula emission and thus we expect higher pulsed 
fraction for pulse profiles derived from XMM-Newton data 
compared to RXTE data. The larger amount of statistical 
uncertainties observed in XMM-Newton data originates 
in lower effective observation time due to the electronic 
chopper set during Timing Mode observation or due to 
low duty cycle in Burst Mode. 

As an example the resulting distribution of the pe- 
riod analysis of the Crab observation of Rev. 0056 is shown 
in Fig. ^. Within statistical errors the pulse periods we 
get from RXTE data agree with the extrapolated radio 
ephemeris marked as solid vertical line. For the XMM- 
Newton data we get a deviation in pulse period of the 
order of \AP/P\ ^ 10"^. This order of magnitude is con- 
sistent for all other objects and readout modes, except 




5.0825»10" 

Trial Period P- 0.033 [s] 

Figures. Resulting -distribution of the period analysis 
of Crab data observed in Timing Mode in Rev. 0056 and 
RXTE PCA. The maximum of each distribution corre- 
sponds to the most significant pulse period. The verti- 
cal line marks the extrapolated radio period, the distribu- 
tion drawn in red corresponds to XMM-Newton Timing 
Mode data. Note that the pulse period derived from XMM- 
Newton data is significantly longer compared to the radio 
period. The width of the distributions can not directly 
be used as error estimate for the pulse period. 



Small Window mode which has a lower time resolution 
(compare Tab. ph. We calculated t^PjP as 



AP 



found 



radi 



radio 



(1) 



where Pradio is the extrapolated radio period. These 
discrepancies are equivalent to a residual velocity com- 
ponent of 0.4-2.2 km/s which might not be taken into 
account during the corrections applied to the photon ar- 
rival times. To cross check our results with EPIC MOS, we 
analyzed data of PSR B1509-58 observed in MOS timing 
mode during Rev. 0137, as well. From this analysis we get 
pulse periods consistent with those derived from EPIC pn 
data. The observed deviations are too large compared to 
the specifications which would imply a AP/P ^ 10~^ and 
prevent a absolute time cross calibration with RXTE for 
the time being. 

To investigate the relative timing in more detail, addi- 
tional calibration observations of the Crab are scheduled. 
These observations will take place at different orbital po- 
sitions of XMM-Newton to be able to separate residual 
orbital effects. 

6. Dead time effects 

In order to quantify the efficiency of the time resolution 
with respect to the dead time for each mode in the EPIC 
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Table 3. Pulse ephemeris resulting from temporal analysis. The expected pulse periods are extrapolated radio periods 
taken from Jodrell Bank monthly pulsar data base and the Princeton Pulsar data base. All radio pulse periods are 
extrapolated periods relative to the epoch given. AP/P is defined as (Pfound — Pradio)/Pradio- 



Object 


RA (J2000.0) 


Dec (J2000.0) 


Mode 


Epoch 


-fradio 


-ffound 


AP/P 










[MJD] 


[ms] 


[ms] 




Crab 


05 34 31.973 


+22 00 52.061 


Timing 


51632.8696 


33.508383 


33.508424 


1.3 X 10-^ 








Burst 


51988.6595 


33.521309 


33.521238 


-2.1 X 10"^ 


B1509-58 


15 13 55.617 


-59 08 08.872 


Timing 


51794.3157 


151.114141 


151.113045 


-7.3 X 10-^ 








Small 


51794.1681 


151.114125 


151.112350 


-1.2 X 10"^ 


B0540-69 


05 40 11.049 


-69 19 55.188 


Timing 


51691.5783 


50.519308 


50.519599 


5.8 X 10"^ 








Small 


51691.0894 


50.519288 


50.519297 


-1.9 X 10"^ 



pn detector we have simulated synthetic time series for 
Timing and Burst Mode. These light curves were analyzed 
using standard time series analysis tools developed at our 
institute. 

The goal of our work was to determine the detection 
probability of a signal in a power spectrum for the fast 
readout modes of XMM-Newton. In all simulations our 
primary goal was to make use of the best possible time 
resolution allowed in the respective mode. Note that some 
of the effects described in the following sections, can par- 
tially be resolved by rebinning the input light curve or the 
power spectrum (PSD), this however in any case degrades 
the time or frequency resolution. A detailed description 
of the simulations and the parameters used is given by 



Kuster et al. (1999) 



6.1. Influence of readout sequence 

For our synthetic light curves we simulated random Pois- 
sonian noise with a time resolution equal to the resolution 
of the observation mode, and with a mean count rate that 
corresponds to the "maximum" count rate allowed for the 
mode according to the telemetry and photon pileup con- 
straints (see Ehle et al. 2001). We subsequently folded the 
data with the readout sequence (dead time function) for 
the respective mode, and computed the power spectrum 
in Leahy normalization. We compared the resulting power 
spectra with those expected for pure white noise. 

The results for Timing and Burst Mode show distorted 
features at large frequencies which are produced by alias- 
ing due to the dead time function of the readout mode (see 
Fig. ^). The first peak appearing in the power spectrum 
corresponds to the life/dead time window in both cases: 
5.91 ms (169.2 Hz) for the Timing Mode and 4.21 ms 
(237.5 Hz) for the Burst Mode respectively. Consecutive 
peaks correspond to multiples of these frequencies. Even 
with rebinning, for these high time resolution modes the 
power spectrum is completely dominated by the window- 
ing above these "critical" frequencies, especially in Burst 
Mode. 

In order to determine whether features in the PSD (e.g. 
QPOs) are detectable above the "critical" frequencies, we 
analyzed the behavior of the power spectrum in the pres- 



ence of a variable source. To represent realistic astronomi- 
cal data, we have simulated a variable X-ray source show- 
ing a quasiperiodic oscillation (QPO) at kHz frequencies 
on top of a red noise spectrum and subsequently folded 
the resulting data with the readout sequence. As the lower 
left PSD panel of Fig. ^ clearly demonstrates, for Timing 
Mode the effect of the sampling rate has no influence on 
the behavior of the QPO in the power spectrum, in the 
sense that the peak in the range of the critical frequen- 
cies can be recognized between the well resolved peaks 
produced by the sampling. On the other hand, for Burst 
Mode, the situation is dramatically different. In this case 
the QPO cannot be detected above the critical frequency 
(see lower right panel of Fig. ^). The sampling rate effect 
completely hides any sign of periodic (or quasiperiodic) 
oscillation in the data. 

6.2. Influence of detector configuration 

In addition to the internal dead time of the pn-CCD due to 
the readout sequence, the configuration of the detector can 
affect the time resolution of an observation as well. In or- 
der to fulfill the telemetry constraints of the EPIC Instru- 
ments especially for bright sources the observer has the 
possibility to use the "electronic chopper" implemented 
in the on board electronics of the EPIC pn camera (see 
[Ehle et al. 2001 ). For a chopper value of N readout 
frames of the CCD are discarded before processing each 
{N + 1) th frame. Thus only the data of each {N + 1) th 
frame is added to the telemetry stream. This adds an ad- 
ditional aliasing window function to the data causing a 
distortion of the PSD. As an example the PSD of an Crab 
observation in Timing Mode is shown in the left panel 
of Fig. ^ During this observation the electronic chopper 
was set to 25, thus reducing the effective life time by the 
same factor. As a result the power spectrum is completely 
dominated by noise above frequencies of 10 Hz. 

Depending on count rate a similar effect is caused by 
counting mode. While the electronic chopper adds a strict 
periodic dead time window, the effect caused by count- 
ing mode generally is quasi-periodic. The observer should 
be aware of these limiting constraints which are of im- 
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Frequency [Hz] 
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Figure 4' Averaged PSDs out of 2500 power spectra assuming Poissonian noise data folded with the readout sequence 
of Timing and Burst Mode (top panel left and right). The same simulations assuming an additional signal of an AR[2] 
process representing a QPO (bottom panel left and right). The mean count rate for the simulations was set to 1500 
cps, with a time resolution of 29.52 jis for Timing Mode and 60000 cps with a time resolution of 72 jis for Burst 
Mode. Note that the QPO can not he detected in Burst Mode. 



portance expecially for bright sources when planning an 
observation. 

7. Conclusions 

The fast modes of the EPIC cameras provide a time res- 
olution which is outstanding for X-ray CCDs. With this 
high time resolution it is possible with X-ray CCD detec- 
tors to resolve pulsed emission even for fast millisecond 
pulsars down to pulse periods of ~ 1 msec. 

Our preliminary analysis of calibration observations of 
young millisecond pulsars clearly demonstrates, that the 
internal timing of the EPIC pn and EPIC MOS camera 
works as expected. However, currently the accuracy of pe- 
riods determined from XMM-Newton observations is only 
of the order of AP/P « 10~^, deviating from an expected 
uncertainty of AP/P ^ 10~^. This prevents a calibration 
of the absolute time for the time being and needs further 
investigation. 



Due to the readout principle of the detector the time 
resolution is limited by aliasing frequencies caused by dead 
time effects. This is of importance especially for bright 
sources and can partially be over come by rebinning the 
light curve or the power spectrum. However, rebinning 
always diminishes the effective time resolution as well. 
In addition the configuration of the detector (electronic 
chopper or too high telemetry rates) can further degrade 
the time resolution by additional periodic (el. chopper) or 
quasi periodic (counting mode) dead time windows. 
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Figure 5. Left: PSD of the Crab observation in Rev. 0056 during which the electronic chopper was set to 25. The PSD 
is totally dominated by noise above 10 Hz (compare Fig. Right: During this pulsar observation the data handling 
unit periodically went to counting mode. The broad peaks below 0.2 Hz are caused by the counting mode. While the 
chopper increases noise at high frequencies, counting mode distorts the PSD at lower frequencies. 



References 

Aranci, G., Cotardo, G., Masaro, D. et al. 1998, XMM-HD- 

LAB-0001, XMM Command and Data Management Unit 

Design Description, ESA 
Briel, U.G., Dennerl, K., Feyberg, M. et al., this volume 
Ehle, M., Breitfellner, M., Dahlem, M. et al. 2001, The XMM 

Users' Handbook, Issue 2.0, ESA 
Kendziorra, E., Bihler, E., Grubmiller, W. et al. 1997, Proc. 

SPIE, 3114, 155 
Kendziorra, E., Colh, M., Kuster, M. et al. 1999, Proc. SPIE, 

3765, 204 

Kirsch, M., Kendziorra, E., Feyberg, M. et al., this volume 
Kuster, M., Benlloch, S., Kendziorra, E. et al. 1999, Proc. 

SPIE, 3765, 673 
Lammers, U. 2001, The XMM ODF Access Layer Library, ESA 
Leahy, D.A., Darbo, W., Eisner, R.F. et al. 1983, ApJ, 266, 

160 

Lyne, A.G., Jordan, C.A., Roberts, M.E. 2001, Jodrell Banks 

Observatory, University of Manchester 
Muno, M.P., Morgan, E.H., Remillard, R.A. et al. 1999, ApJ, 

527, 321 

Taylor, J.H., Manchester, R.N., Lyne, A.G. 1993, ApJS, 88, 
529 



